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Contract/Grant  Title:  Quantum  Transport  and  Optoelectronics  in 
Gapped  Graphene  Nanodevices 

Contract/Grant  #:  FA9550-1 1-1-0225 

Final  Report 

During  the  period  funded  by  this  grant,  we  investigated  the  fundamental  physical  processes  of 
light-matter  interactions  in  van  der  Waals  crystals,  including  graphene,  boron  nitride,  transition 
metal  dichalcogenides  and  heterostructures  formed  from  them.  Our  research  leads  to  deep 
understanding  of  novel  physics  embedded  in  the  low  dimensionality,  such  as  ultrafast  electron- 
electron  interaction,  unconventional  electron-phonon  coupling  and  strong  light  confinement. 
This  understanding  paves  the  way  for  creating  high-perfonnance  optoelectronic  and  energy 
harvesting  devices  with  state-of-the-art  functionalities,  equipped  by  these  innovative  materials 
and  unique  device  architectures.  This  grant  has  resulted  in  12  high  profde  publications,  including 
two  papers  in  Science^  1,  2],  two  in  Nature  Physics[ 3,  4],  two  in  Nature  Nanotechnology[5,  6], 
two  in  Nature  Communications^! ,  8],  one  in  Physical  Review  Letter s\ 9],  two  in  Nano  Letter s[  10, 
11]  and  one  in  Advance  Materials[\2\.  Below  we  highlight  a  few  major  works  funded  by 
AFOSR. 

1.  Hot  Carrier-Assisted  Photoresponse  in  Intrinsic  Graphene[l,  6,  9] 

The  photoresponse  of  materials,  which  determines  the  performance  of  optoelectronic  devices,  is 
governed  by  energy  relaxation  pathways  of  photoexcited  electron-hole  pairs:  energy  transferred 
to  the  lattice  is  lost  as  heat,  while  energy  retained  in  the  electronic  subsystem  can  be  used  to 
drive  an  optoelectronic  circuit.  In  graphene,  energy  relaxation  pathways  are  strongly  affected  by 
a  vanishing  electronic  density  of  states,  which  creates  a  bottleneck  that  limits  energy  transfer  into 
the  lattice.  As  a  result,  the  photogenerated  carrier  population  remains  hot  while  the  lattice  stays 
cool. 

In  graphene,  hot  carriers  should  play  a  key  role  in  the  optoelectronic  response.  Initial  studies 
observed  photocurrent  response  both  at  graphene-metal  contacts  and  monolayer-bilayer  graphene 
junctions,  yet  the  underlying  mechanism  was  under  debate  between  the  photovoltaic  (PV)  effect, 
where  the  photogenerated  electron-hole  pairs  are  separated  by  a  built-in  electric  field  [13-15], 
and  the  photothennoelectric  (PTE)  effect,  driven  by  electron  temperature  gradients  [16],  In  our 
work,  we  carried  out  scanning  photocurrent  measurements  of  highly  controllable  dual-gated 
graphene  p-n  junction  devices  (Fig.  3a)  to  resolve  this  debate.  Tuning  the  bottom-  and  top-gate 
voltages,  Vbg  and  Vtg  respectively,  allowed  independent  control  of  carrier  density  of  electrons 
and  holes  in  each  region  to  form  a  p-n  junction  interface  at  the  middle  of  the  device  (Fig.  3b) 
[17-19].  Figure  3d  shows  a  photocurrent  image  obtained  by  scanning  a  focused  (1  pm)  laser  spot 
over  the  device  (Fig.  3c)  while  measuring  the  current  I. 

Our  measurements  indicated  that  hot  electronic  carriers  dominate  the  intrinsic  optoelectronic 
response  of  graphene  in  the  linear  optical  power  regime.  The  hot  carrier  regime  manifests  as  a 
strong  PTE  effect  that  results  in  a  striking  six-fold  sign-alternating  photocurrent  pattern  as  a 
function  of  gate  voltages  (Fig.  3e),  qualitatively  different  from  the  PV  effect  which  would  only 
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have  a  two-fold  sign-alternating  photocurrent  pattern.  Additionally,  the  spatial  pattern  and  the 
charge  density  dependence  of  the  photoresponse  established  a  strong  connection  between 
thermal  energy  transport  and  electronic  charge  transport.  In  line  with  this  finding,  the  graphene 
p-n  junction  device  can  be  operated  as  a  novel  photo-thennometer  that  measures  local  electron 
temperature.  In  collaboration  with  the  Koppens  group  in  ICFO,  we  carried  out  a  time-resolved 
photovoltage  measurement  and  resolved  a  voltage  generation  time  of  less  than  50  femtoseconds. 
This  ultrafast  response  was  further  demonstrated  by  electrically  measuring  the  pulse  duration  of 
a  sub-50  femtosecond  laser  pulse.  The  above  results  have  been  published  in  Science  [  I  ]and 
Nature  Nanotechnology  [6]. 

Continuing  this  line  of  work,  we  performed  a  systematic  temperature  dependence  of  the  PTE 
current  to  detennine  the  dominant  hot  carrier  cooling  channels  in  graphene.  A  persisting  six-fold 
pattern  at  all  temperatures  (4  -  300  K)  indicates  that  the  PTE  effect  always  governs  the 
photocurrent  response,  allowing  us  to  track  the  induced  electron  temperature  and  probe  hot 
electron  cooling  mechanisms.  We  observed  a  pronounced  non-monotonic  temperature 
dependence  of  the  PTE  photocurrent  (Fig.  3f),  which  we  understand  as  the  competition  between 
two  hot  electron  cooling  pathways,  i.e.,  momentum-conserving  normal  collisions  that  dominate 
at  low  temperatures  and  disorder-assisted  supercollisions  that  dominate  at  high  temperatures  [20- 
22].  The  peak  temperature,  where  the  total  cooling  rate  is  minimal,  strongly  depends  on  charge 
carrier  density  and  disorder  strength,  thus  allowing  for  an  unprecedented  way  of  controlling 
graphene's  photoresponse.  This  work  was  selected  as  an  Editors  ’  Suggestion  in  Physical  Review 
Letters  [9]. 


Figure  1.  Hot  carrier-assisted  photoresponse  in  graphene  and  its  dependence  on  sample 
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temperatures,  a,  Photocurrent  measurement  scheme  with  scanning  laser  excitation,  b, 
Schematic  of  band  alignment  for  a  monolayer  graphene  p-n  junction,  c,  Optical  micrograph  of  a 
dual-gated  graphene  device  incorporating  BN  as  the  top  gate  dielectric,  d,  Spatially  resolved 
photocurrent  map  with  laser  wavelength  850  mn.  White  lines  mark  the  location  of  gold  contacts 
and  top  gate  electrode,  e,  Photocurrent  of  the  graphene  p-n  junction  versus  back-gate  and  top- 
gate  voltages,  measured  with  the  laser  fixed  at  the  p-n  interface  (black  triangle  in  d).  The  six¬ 
fold  pattern  of  the  photovoltage  sign  changes  provides  evidence  that  the  PTE  effect  rather  than 
the  PV  effect  dominates  the  intrinsic  photoresponse  in  graphene,  f,  The  photovoltage  as  a 
function  of  the  sample  temperature  (nonnalized  to  the  maximum  value)  at  particular  gate  voltage 
values  (black  circles  in  e)  along  the  white  dashed  line,  where  the  chemical  potentials  are  equal 
but  opposite  in  signs  at  the  two  parts. 

2.  Tuning  Electron  Thermalization  Pathways  in  a  van  der  Waals 
Heterostructure  [3] 

Immediately  after  photoexcitation  of  an  optoelectronic  device,  energetic  electrons  scatter  with 
other  high-energy  and  ambient  charge  carriers  to  form  a  thennalized  hot  electron  gas,  which 
further  cools  by  dissipating  excess  energy  to  the  lattice.  Due  to  the  short  distance  travelled  by 
charge  carriers  between  electron-electron  scattering  events  in  solids  [23],  equilibration  among 
the  electrons  occurs  on  the  tens  of  femtoseconds  to  picosecond  time  scales  [24,  25].  In  graphene, 
a  low-dimensional  material  with  much  enhanced  Coulomb  interaction[26],  electron 
thennalization  is  known  to  occur  on  extremely  fast  time  scales  (<30  fs)  [6,  27-29],  reflecting  the 
extremely  short  transit  length  between  scattering  events.  Most  analyses  of  graphene  have, 
therefore,  treated  its  electrons  to  be  instantaneously  thennalized  [1,9,  20,  21],  and  slightly  non- 
thennal  electronic  behavior  has  thus  far  only  been  reported  in  pump-probe  experiments  with 
ultrashort  (-10  fs)  laser  pulses  and  low  excitation  density  [30,  31].  Due  to  such  short  time 
(femtosecond)  and  length  (nanometer)  scales,  it  is  challenging  to  detect  and  control  the 
thermalization  process  in  graphene,  or  more  generally,  in  any  solid-state  systems. 

We  take  a  novel  approach  to  probe  and  manipulate  electron  thermalization  in  graphene  by 
introducing  a  new  energy  transport  channel  that  competes  with  the  thermalization  process.  Such 
an  additional  dynamical  pathway  is  realized  in  a  van  der  Waals  heterostructure  that  consists  of  a 
graphene-BN-graphene  stack  (Fig.  4a-c).  In  this  layered  structure,  the  photoexcited  electrons  in 
one  graphene  layer  can  travel  vertically  to  the  other  graphene  layer  through  the  very  thin  middle 
BN  layer  (blue  dashed  arrow  in  Fig.  4b).  Given  the  close  proximity  of  these  layers,  interlayer 
charge  transport  can  occur  on  extremely  fast  time  scales,  competing  directly  with  the  intralayer 
thermalization  process  (red  arrows  in  Fig.  4b). 

In  our  experiment,  we  have  observed  such  a  competition  by  measuring  the  interlayer 
photocurrent  under  different  bias  and  excitation  conditions  (color  map  shown  in  Fig.  4d-e).  In 
particular,  we  observe  distinct  power-dependences  of  the  photocurrent  in  the  different  regimes, 
as  shown  in  Fig.  4f.  Region  A  (low  photon  energy  and  low  bias  voltage)  shows  highly- 
superlinear  power  dependence,  whereas  Region  B  (high  photon  energy  and  high  bias  voltage) 
shows  linear  power  dependence.  These  two  distinct  power  dependences  are  attributed  to  two 
distinct  photocurrent  generation  mechanisms  in  our  structure:  thermionic  emission  and  direct 
tunneling.  By  adjusting  the  interlayer  bias  voltage  and/or  the  excitation  photon  energy,  we  can 
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tune  the  interlayer  charge  transport  to  be  slower  or  faster  than  the  intralayer  thermalization,  thus 
tuning  the  overall  thermalization  process. 

Our  experiments  not  only  provide  valuable  insight  into  the  electron  dynamics  of  graphene,  but 
also  demonstrate  a  new  means  to  manipulate  electron  thermalization  in  low-dimensional 
materials.  With  appropriate  modeling,  we  have  deduced  the  thermalization  time  in  graphene  (~10 
fs)  from  a  time-averaged  transport  experiment  (Fig.  4f  black  dashed  lines).  Similar  experimental 
methods  can  readily  be  extended  to  other  van  der  Waals  heterostructures  to  probe  their  ultrafast 
electronic  processes.  More  generally,  given  the  central  role  that  electron  thermalization  plays  in 
energy  transport,  the  tunability  on  femtosecond  time  scales  demonstrated  in  our  research  opens 
up  a  new  range  of  in-situ  functionality  for  novel  optoelectronic  devices.  This  work  was  recently 
published  in  Nature  Physics[3]. 
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Figure  2.  Interlayer  photocurrent  of  a  G-BN-G  device,  a,  Schematic  of  a  G-BN-G  device 
under  optical  excitation,  b,  Schematic  of  intralayer  thermalization  and  interlayer  transport  of  the 
optically  excited  carriers,  c,  Optical  image  of  a  G-BN-G  device,  which  consists  of  a  top 
exfoliated  graphene  layer  (white  line),  a  14-nm-thick  BN  flake,  and  a  bottom  graphene  layer 
grown  by  chemical  vapor  deposition  (dashed  red  line),  d-e,  Scanning  images  of  interlayer 
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photocurrent  at  interlayer  bias  Vt  =  -0.5  and  0.5  V.  f,  A  color  map  of  y,  the  exponent  of  the 
power-dependence  of  the  photocurrent,  as  a  function  of  Vb  and  tieo.  The  color  scale  is 
customized  to  make  the  area  with  y  >  1  (Region  A)  and  y  =  1  ±  0.02  (Region  B)  appear  red- 
yellow  and  blue,  respectively.  The  black  dashed  lines  are  theoretical  contours  corresponding  to 
tunneling  time  of  2,  7,  100  and  1000  fs,  predicted  by  our  model  described  in  Ref.  [3].  g, 
Schematics  depicting  the  thennionic  emission  and  the  direct  carrier  tunneling  as  the  dominant 
photocurrent  mechanism  for  Region  A  and  B  in  (f),  respectively.  All  measurements  were  carried 
out  with  a  supercontinuum  laser  at  T=  100  K. 


3.  Collective  Excitations  in  Hexagonal  Boron  Nitride  and 
Graphene  [2,  5,  7] 

Hexagonal  boron  nitride  has  been  widely  used  as  high  quality  substrate  for  low-dimensional 
materials  and  also  as  a  capping  layer  to  protect  these  materials  from  atmospheric  contaminants 
and  oxidation[32,  33],  In  addition,  heterostructures  of  BN  with  other  2D  materials  can  form 
superlattices  (or  moire  patterns),  which  strongly  alter  the  combined  electronic  properties[34-36]. 
As  a  polar  material,  light  can  strongly  couple  to  BN’s  optical  phonon  resonance,  and  the  layered 
nature  can  lead  to  novel  properties  of  this  resonance.  We  collaborated  with  Prof.  Dimitri  Basov’s 
group  at  UC  San  Diego  to  investigate  phonon-polariton  propagation  in  this  layered  material 
using  an  infrared  nanoimaging  technique  in  a  scattering-type  near  field  scanning  optical 
microscope  (s-SNOM). 

We  launched,  detected,  and  imaged  polaritonic  waves  in  real  space  and  altered  their  wavelength 
by  varying  the  number  of  crystal  layers  (Fig.  5a-g).  The  measured  dispersion  of  polaritonic 
waves  was  shown  to  be  governed  by  the  crystal  thickness  according  to  a  scaling  law  that  persists 
down  to  a  few  atomic  layers.  We  also  discovered  that,  at  certain  frequencies,  BN  exhibits 
hyperbolic  properties:  the  axial  and  tangential  permittivities  have  opposite  signs.  In  such 
circumstances  light  propagation  is  unusual,  leading  to  novel  optical  phenomena.  We  reported 
infrared  nano-imaging  experiments  demonstrating  that  BN  can  act  as  a  “hyper-focusing  lens” 
and  as  a  multi-mode  waveguide  in  mid-infrared  frequencies.  The  lensing  is  manifested  by 
subdiffractional  focusing  of  phonon-polaritons  launched  by  pre-made  metallic  patterns 
underneath  the  BN  crystal  (Fig.  5h-i).  Our  work  opens  new  opportunities  for  anisotropic  layered 
insulators  in  infrared  nanophotonics  complementing  and  potentially  surpassing  concurrent 
artificial  hyperbolic  materials  with  lower  losses  and  higher  optical  localization.  This  work  was 
published  in  Science  [2]  and  Nature  Communications  [7]. 

In  addition,  we  have  implemented  and  investigated  the  tunable  hyperbolic  response  in 
heterostructures  comprised  of  a  monolayer  graphene  deposited  on  hexagonal  boron  nitride  (G- 
BN)  slabs.  Electrostatic  gating  of  the  graphene  layer  enables  electronic  tuning  of  the  phonon- 
polariton  properties  in  BN.  The  tunability  originates  from  the  hybridization  of  surface  plasmon- 
polaritons  in  graphene  to  hyperbolic  phonon-polaritons  in  BN,  which  we  examined  via  nano-IR 
imaging  and  spectroscopy.  These  hybrid  polaritons  possess  a  combination  of  properties  from 
plasmons  in  graphene  and  phonon-polaritons  in  BN.  Therefore,  G-BN  structures  fulfill  the 
definition  of  an  electromagnetic  metamaterial,  since  the  combined  properties  of  these  devices  is 
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not  revealed  by  the  constituent  elements  alone.  Our  results  uncover  a  practical  approach  for 
realization  of  nano-photonic  metamaterials  by  exploiting  the  interaction  of  distinct  types  of 
polaritonic  modes  hosted  by  different  constituent  layers  in  van  der  Waals  heterostructures.  This 
work  was  published  in  Nature  Nanotechnology  [5]. 
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Figure  3.  Real-space  imaging  of  surface  phonon-polaritons  on  BN  and  a  BN-based 
‘hyperlens’.  a,  Schematics.  Arrows  denote  the  incident  and  back-scattered  IR  light.  Concentric 
yellow  circles  illustrate  the  phonon  polariton  waves  launched  by  the  AFM  tip  and  reflected  by 
the  two  edges  of  a  tapered  BN  crystal,  (b  and  d  to  f)  IR  near-field  images  of  the  normalized 
amplitude  signal  and  taken  at  different  IR  frequencies  (BN  thickness  in  b  to  f  d  =  256  nm).  c, 
Simulation  of  the  phonon-polariton  interference  pattern,  g,  Phonon-polaritons  probed  in  three- 
layer  (left)  and  four-layer  (right)  BN  crystals.  The  white  dashed  line  tracks  the  BN  edges 
according  to  the  AFM  topography.  Scale  bars  indicate  800  nm.  h-i,  Use  BN  as  a  hyper-lensing 
material  to  image  the  sub-diffractional  structure  (letter  ‘UCSD’  and  ‘MIT’)  underneath. 
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